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ABSTRACT 


We  have  found  a  method  for  obtaining  algebraic  approximations  to  certain  integrals 
that  appears  to  be  independent  of  known  methods,  such  as  Taylor,  orthonormal  functions, 
Laplace’s  method,  etc.  This  method  has  the  advantage  of  sometimes  producing  much  more 
compact  forms  of  the  integral.  It  often  gives  results  where  the  others  may  fail  and  in  many 
cases  the  solutions  are  globally  valid  over  the  full  range  of  each  parameter.  It  sometimes 
allows  for  a  trade-off  between  the  accuracy  and  the  functional  complexity  of  the  solution. '' 


Specific  problems  in  the  estimation  of  extinction  from  obscurants  and  other  aerosols 
are  for  the  first  time  completely  resolved  analytically.  These  analytical  solutions  are  ap¬ 
plicable  to  infrared  and  millimeter  wave  smokes  as  well  as  biological  and  natural  aerosols.  Ct 


RESUME 

Nous  avons  trouve  une  methode  qui  permet  d’obteuir  des  approximations 
alg^briques  a  certaines  int^ales.  Cette  methode  nous  apparait  ind^pendante  des  autres 
m^thodes  connues  comme  celle  de  Taylor,  de  Laplace,  ou  celle  des  expansions  en  fonctions 
orthonormales.  Dans  certains  cas,  cette  methode  a  I’avantage  de  produire  des  formes  beau- 
coup  plus  compactes  de  Tmt^ale,  de  donner  an  r^ultat  alors  que  les  autres  m^thodes 
4chouent  et  de  produire  dect  solutions  globalement  valides  sur  toute  l’4chelle  des  valeurs  de 
chaque  parametre.  Elle  permet  souvent  de  balancer  la  pr^sion  du  r^ultat  et  la  complexity 
fonctionnelle  de  la  solution. 

Certains  probiymes  qui  concement  I’extinction  causye  par  les  agents  obscurcissants 
et  par  les  ayrosols  ont  yty  compiytement  solutionnys  analytiquement  pour  la  premiyre  fois. 
Ces  solutions  analytiqnes  sont  applicables  aux  fumyes  actives  dans  les  rygions  infrarouges 
et  millimytriques  ainsi  qu’aux  ayrosols  d’oHgine  biologique  et  naturelle. 
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Any  complex  military  system  or  scenario  that  involves  electro-optics  requires  the 
solution  to  one  or  many  diffraction  problems.  This  includes  propagation  of  radiation  through 
natural  aerosols,  smokes /obscurants,  design  and  evaluation  of  high  precision  optical  systems, 
novel  laser  cavities,  etc. 

Diffraction  problems  nearly  always  involve  the  evaluation  of  a  certain  class  of  in¬ 
tegrals.  These  integrals  are  usually  intractable  and  thus  must  be  computed  numerically. 
Hence,  this  can  be  not  only  computationally  expensive  but  there  is  in  general,  significant 
loss  of  physical  understanding. 

This  report  describes  a  technique  that,  for  many  diffraction-type  integrals,  gives  an 
approximate  but  efficient  analytic  solution.  This  reduces  the  computational  effort  consid¬ 
erably  while  retaining  the  physical  insight  essential  for  system  design  or  scenario  analysis. 

Specific  integrals  involved  in  the  estimation  of  extinction  from  obscurants  and  other 
aerosols  are  for  the  first  time  completely  resolved  analytically.  These  analytical  solutions 
are  applicable  to  infrared  and  millimeter  wave  smokes  as  well  as  to  biological  and  natural 
aerosols. 

Although  all  examples  given  in  this  work  originate  from  the  diffraction  of  electro¬ 
magnetic  radiation  by  aerosols,  the  integration  technique  is  not  necessarily  restricted  to  this 
area  of  physics.  We  therefore  expect  the  technique  to  be  useful  in  other  fields.  It  is  at  this 
stage  a  possible  solution  looking  for  other  promising  problems. 
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There  are  some  areas  in  physics  where  an  analytic  approximation  to  an  integral  is 
highly  desirable  even  though  it  is  inexact.  This  is  particulariy  true  when  the  physical  theory 
is  known  to  be  only  an  approximation  itself.  A  striking  example  are  integrals  that  occur 
in  the  Fresnel-Huygens  or  Kirchhoif  form  of  scalar  diffraction  theory  which  is  applicable 
to  both  optics  and  acoustics.  Aerosol  scattering  problems  provide  another  example.  One 
wishes  for  an  approximate  analytic  solution  in  order  to  retain  at  least  one  of  the  follov/ing: 
physical  insight,  correct  asymptotic  behaviour,  computatknal  speed  and/or  computational 
robustne».'<. 

We  have  found  a  method  for  obtaining  algebraic  ^proximations  to  certain  integrals 
that  appears  to  be  independent  of  known  methods,  such  as  Taylor,  orthonormal  functions, 
method  of  stationary  phases,  steepest  descent,  etc.  This  method  has  the  advantage  of 
often  producing  much  more  compact  forms  of  the  integral  and  sometimes  producing  results 
where  the  others  may  fail.  It  sometimes  allows  for  a  trade-off  between  the  accuracy  and 
the  functionrd  complexity  of  the  solution.  For  the  practical  purposes  of  this  report,  when 
we  refer  tc  an  analytic  solution  we  mean  that  the  result  integration  must  be  represented 
in  terms  of  combinations  of  the  elementary  and  special  functions  of  mathematical  physics. 

The  integrals  that  can  be  approximated  by  this  method  are  the  same  as  those  found 
in  integral  tables  but  the  arguments  of  some  of  the  functions  in  the  kernel  are  much  more 
complicated.  Ai  will  be  shown,  the  main  limitation  of  the  technique  is  that  there  is  no  a 
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priori  error  estimation  nor  is  the  method  as  yet  fully  automated. 

This  work  was  performed  at  DREV  between  July  and  December  1990 ’under  PSC 
32A,  EO/IR  Protection  of  Land  Vehicles. 


2.0  METHOD 


This  chapter  gives  the  general  mathematical  framework  required  to  implement  the 
integration  technique. 


2.1  Outline  of  Concept 


Suppose  one  has  the  following  integral  to  perform: 

/  =  [1] 

where  g,  h  and  /  are  any  continuously  differentiable  functions  and  rjj  and  a/,  ^  param¬ 
eters  independent  of  z  whose  domain  can  cover,  in  general,  the  complex  plane.  In  order  to 
refer  to  the  whole  set  of  parameters  we  will  define  p  =  {A,,r7j,<7ik}>  As  it  stands,  [1]  is  a 
genera’,  integral  with  a  function  that  has  a  potentially  very  complicated  argument  h{tjj,x). 
Most  integrals  that  have  analytic  solutions  have  rather  simple  arguments.  The  following 
method  extends  the  list  of  analytically  soluble  integrals  to  a  list  that  includes  analytically 
approximated  integrals. 
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We  start  oif  by  assuming  that  a  function  u(qi,h)  can  be  found  such  that 

ri^(n„y)  pi  ,  , 

/  /  f{ffk,x)dx  [2] 

is  an  acceptable  approximation  valid  over  the  range  of  integration.  Here  y  ranges  between 
a  and  b,  and  the  9/  are  n  fitting  parameters.  The  number  of  fitting  parameters  n  will 
generally  depend  on  the  functional  form  chosen  for  u.  We  next  apply  Leibniz’s  theorem  for 
differentiation  of  an  integral  on  [2]  to  get: 


d  /*(»»> .»)  d  /» 

—  /  u(qi,h)dh  as  -r-  /  /(fft,i)di  which  gives 

ay  • 


a(qt,h(y))dh  «  /(<r*,»)dy. 


[3] 


Multiplying  each  side  of  [3]  by  g(h)  and  integrating  we  get: 

/«  g{h)u{h)dh  [4] 

where  a  =  h{r}j,a)  and  ^  =  h{qj,b)  to  simplify  the  notation.  Whether  [2]  is  an  acceptable 
approximation  will  depend  on  the  desired  accuracy  of  [4]. 


It  IS  required  that  [4]  be  analytically  integrable.  It  is  often  preferable,  but  not 
necessary,  that  both  the  integrals  in  [2]  also  be  analytically  integrable.  It  is  clear  that  when 
the  approximation  [2]  is  exact  this  procedure  is  just  a  variable  substitution. 


Equation  [2]  can  be  further  generalized  to 

rO. 


[5] 
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with  weighting  constants,  keeping  in  mind  that  for  each  k 

gik)uuih)dh  [6] 

JQ 

must  be  analytically  integrable.  Thus  for  each  function,  u*,  that  satisfies  the  integrability 
conditions  in  [6],  an  additional  functional  form  is  found  that  may  improve  the  approximation 
in  [5]. 

The  concept  above  will  be  made  clearer  by  an  outline  of  some  possible  procedures 
and  several  examples  taken  from  actual  physical  problems. 

2.2  Outline  of  a  Procedure 

In  the  following  subsections  we  will  outline  the  steps  of  one  proced’:.re  for  putting 
into  practice  the  previous  concept. 

2.2.1  Finding  Candidate  u 

In  deriving  a  practical  procedure  for  implementing  the  above  method,  it  is  useful  to 
prune  the  set  of  psesible  forma  of  u.  The  strongest  constraint  is  the  integrability  condition 
[4]  or  [5].  One  way  of  determining  part  of  the  allowed  set  of  functional  forms  for  u  is  by 
looking  up  a  table  cX  definite  integrals  that  contain  products  of  g  and  other  functions.  A 
more  general  way  is  to  consider  the  small  list  of  general  integrals  (Refs.  1-4).  This  list 
can  be  used  to  rapidly  survey  almost  all  known  analytic  definite  integrals.  Although  these 


f 
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integrals  axe  complicated,  they,  at  the  very  least,  serve  as  indicators  as  to  the  existence  of 
possible  forms  of  u.  Examples  from  this  list,  two  of  which  we  will  use,  are  given  below. 


••  • » 

?(./  +  «)  i.,  h,. 

•  •  • » 

•  «  •  f 

Op,  hA 

ft**  j»  ) 

[7] 

^  i~'(z  -  l)'jFi(«  +  <7-p,A  +  <r-p;(7;l  ^ 

" 

. . . ,  fly  1 

* " »  ft*  / 

1  dx 

-  A,  i, .  6,j 

[8] 

_  r(<y)  r(/>)  f  (o). 

^^C+m+n:DiD>  \  (^e). 


(6);  (&') 
(d);  i<f) 
2. 


-*r) 


dx 


r{<T  +  0) 


m-f  n 

(ft);  (ft') 
(‘^);  {<f) 


m'"a 


»• 

.  •  • , 

n*/3 


fn<f  n 


(m  +  ’ (m  +  n)"‘+’* 


)• 


[9] 


Equations  [7]  and  [9]  are  examples  of  general  Euler  type  integrals.  In  [7],  jFj  is  the 
Gauss  hypergeometric  function  and  is  Meijer’s  G-function.  In  [9]  is  a  Kamp4 

de  F^riet  function.  This  double  variable  hypergeometric  function  and  its  properties  are 
discussed  at  length  in  (Ref.  5)  and  more  recently  in  (Ref.  4).  One  very  important  property 
of  the  Kamp4  de  F^riet  functions,  that  we  will  use  iu  an  example,  occurs  when  C  and  A  are 
both  zero  in  which  case  they  become  just  the  product  of  two  generalized  hypergeometric 
functions  of  a  single  variable: 


gAB,B>  ((a): 


(ft); 

(dy, 


(ft') 

id') 


((j)  I') 


Iio| 
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Equation  [8]  is  a  general  Weyl  integral  over  the  product  of  a  Gauss  hypergeometric 
function  and  a  Meijer’s  G-function.  The  detailed  conditions  of  validity  of  these  integrals 
will  not  be  stated  here  so  as  not  to  sidetrack  the  main  issues.  They  can  however  be  found 
in  the  cited  references. 

With  these  integrals  the  functions  g  and  u,  plus  the  appropriate  integral  transform 
can  be  identified.  Hence  j  or  u  could  be  either  a  Meijer’s  G-function,  a  generalized  hyper- 
geometric  function,  simple  functions  (such  as  i^,  (1  -  i)**,  c*  or  sin(i)),  etc.  Thus  most 
functions  occun-ing  in  mathematical  physics  are  covered. 

Additional  qualitative  considerations  of  the  functional  form  of  the  right-hand  side 
of  [2]  can  further  reduce  the  set  of  candidate  forms  for  «.  Some  of  the  more  obvious  ones 
are  the  monotonicity,  curvature  and  asymptotic  behavior  ir  both  large  and  small  limits.  A 
more  subtle  condition  is  the  weighting  provided  by  g  in  (4). 

2.2.2  Obtaining  the  Fitting  Parameters  qi 

The  fitting  parameters,  (in  general,  but  not  always,  functions  of  p,  a  and  6),  in 
[2]  are  obtained  by  some  appropriate  minimizing  procedure  that  reduces  the  error  in  the 
approximation  [4].  This  can  be  stated  in  one  form  as: 

Min  irif.  '  ^  9{n)u{qi,h)dh-  j(h(z))/(i)  deej  dy 


(11] 
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where  the  Min  is  understood  to  be  perfonned  for  a  particular  set  of  p.  The  process  is 
repeated  over  the  parametric  space  of  interest  in  p.  If  the  level  of  error  is  not  satisfactory, 
another  candidate  form  for  u  must  be  tried  either  in  [2]  or  [5].  In  the  latter  case,  we  note 
that  fitting  of  ait  is  also  required. 

Possible  minimizing  strategies  could  Include  numerical  methods  such  as  quasi- 
Newton,  conjugate  gradient  or  simplex,  etc.  These  methods  initiady  do  not  produce  analytic 
expressions  of  qi  in  terms  of  />,  p,  a  and  b.  However,  provided  that  the  dependency  i>  not 
too  complicated,  acceptable  functional  approximations  may  be  found. 

The  above  process  is  ccmplicated,  can  be  computationally  demanding  and  ulti¬ 
mately  could  be  intractable.  However,  th:  examples  in  the  next  chapter  show  that  often 
the  situation  can  be  considerably  simplified  by  relaxing  the  constraint  equation  [11].  Indeed 
it  will  be  shown  that  even  a  trial  and  error  method  for  finding  «  can  often  suffice.  The 
constraint  condition  [11]  is  also  reduced  by  first  assuming  that  ;  is  a  constant,  and  then 
fitting  at  only  a  small  set  of  points,  yi,  equal  to  the  number  of  fitting  parameters  n.  Thus 
[11]  becomes 

u(qi,h)dh-  I  /{x)dx\  =  0,  [12] 

The  function  g  is  only  peripherally  considered  when  choosing  yi  and  u.  As  will  be  shown 
by  examples,  [12]  sometimes  can  be  solved  analytically,  i.e.  qi  in  terms  of  h,  p,  a  and  6. 
This  simplified  procedure  has  been  shown  to  work  for  functions  with  simple  forms  such  as 
monotonic  ones,  whose  curvature  is  easily  modelled. 
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3.0  EXAMPLES 

In  this  chapter  we  will  give  several  examples  of  the  application  of  the  method  to 
various  integrals  that  occur  in  optical  scattering  theory.  None  of  these  can  be  transformed 
into  integrals  that  have  closed  form  solutions  in  terms  of  special  functions.  Moreover,  since 
the  parameters  in  each  kernel  can  arbitrarily  take  on  values  from  0  to  oo,  Taylor  series 
expansions  and  asymptotic  methods  such  as  Laplace’s  or  steepest  descents  will  not  work 
because  of  convergence  problems.  Since  part  of  the  kernel  has  a  complicated  argument, 
any  expansion  in  sets  of  orthogonal  functions  will  also  involve  this  argument  and  hence  will 
most  likely  make  any  possible  solution  undesirable.  In  all  three  cases,  the  solution  obtained 
is  valid  globally  over  the  complete  range  of  each  parameter  in  p.  This  last  point  will  be 
illustrated  by  extensive  tables. 

3.1  Example  1 

The  following  example  occurs  in  the  study  of  scalar  scattering  from  randomly  ori¬ 
ented  inhnite  aspect  discs. 

Qi  =  ::(l-Re(/i))  with 

Id  =  [13] 

where  u;(i)  =  2j3[(m*  -  -  (1  —  s  is  the  particle  size  parameter,  a  real  number 

varying  between  0  and  oo,  and  m  is  the  refractive  index  and  is  hence  complex  with  Re(Tn) 
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varying  between  1  and  oo  and  Im(m)  varying  between  0  and  — oo.  Hence  the  magnitude  of 
a;(x)  varies  from  0  to  oo  as  well. 

Following  the  notation  of  the  previous  chapter,  h()7i,r;a,x)  =  w,  g{h)  =  e~^,  /(x)  = 
1,  a  =  0,  and  4=1.  The  parameters  are  obviously  rji  =  s  and  %  =  m  and  thus  p  = 

The  functional  form  of  u  is  constrained  by  the  integrability  condition  [4].  As  g{h)  = 
e~*,  the  simplest  nontrivial  possible  forms  of  u{h)  are  limited  to  ce~^^  or  1/(4+  where  c 
and  d  are  constants.  The  term  “simplest”  is  defined  here  by  the  order  of  the  hypergeometric 
function  used.  For  these  forms  e*  is  a  qFq  and  1/(4  +  c)  is  a  iFo,  all  other  hypergeometric 
functions  are  much  more  complicated  in  their  behavior,  computation  and  in  the  integrals 
that  would  result  from  this  procedure.  Of  the  two  forms  identified  above,  ce”***  gives  rise 
to  a  set  of  transcendental  equations  that  needs  to  be  solved  in  order  to  obtain  c  and  d.  We 
thus  choose  the  form  1/(4  +  c)‘^  which  produces  readily  solvable  equations.  This  form  gives 
the  approximation  of  I  in  terms  of  exponential  integrals  of  order  d. 


Let 


rvi 

Jiyi)=  /  . 

Jo 


/(x)  dx  =  gi. 


[14] 


Since,  over  the  range  of  physical  interest  Re(m)  >  1,  4  is  a  simply  behaved  function 
of  yi,  we  can  find,  in  this  case,  a  transformation  relating  4  to  J.  Here,  this  transformation 
u  can  be  of  the  form  stated  above.  Thus  we  have 


1 

(4  +  c)K-‘)- 


[15] 


UNCLASSIFIED 

10 


We  have  empirically  found  that  d  =  2  is  a  reasonably  good  approximation  to  J. 
Furthermore  since  the  above  form  allows  a  maximum  of  three  fitting  parameters,  n  =  3. 

Thus  we  find 

J(yi)  =  yi«^ ..  1  =  1,2,3.  (16] 

Kvi)  +  ft 

This  set  of  equations  can  be  expressed  as  three  linear  equations  which  are  readily  solved 
for  gi.  The  values  of  y/  were  chosen  so  that  the  equations  are  simple,  i.e.  yi  -  0,  yj  =  1/2 
and  y3  =  1.  Other  sets  of  yj  could  be  used  achieving  a  slightly  better  approximation  at  the 
expense  of  algebraic  complexity.  With  the  chosen  set  of  yj,  the  solution  to  [16]  is 


ft  =  +  ft)ft 

h(0)h(l/2)  -  2A(0)/i(l)  +  h(l/2}h(l) 
^  h(0)-2h(l/2)  +  h(l) 

_  /t(l)-b  ft 

^  "  Ml) -MO)' 


[17] 


In  the  above  we  have  simplified  the  notation  for  hCft.ft.y/)  to  h{yi)  for  clarity. 


With  the  above  equations  we  can  substitute  into  [4]  to  get: 

/■Ml)  e-A 

U  »  /  77— — ^dh 

Jh(o)  (*  +  ft)^ 

=  ,  ■^[g^(Ml)  +  ft)  E^{h{0)  +  q2)] 

L  Mi)  +  ft  MO)  +  ft  J 


[18] 


Note  that  Ej  is  the  second  order  exponential  integral.  Thus  the  integral  [13]  has 


been  approximated  as  the  difference  of  two  second  order  exponential  integrals. 
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Table  I  presents  a  comparison  between  the  exact  numerical  calculation  of  and 
the  approximation  [18]. 

As  seen  in  Table  I,  when  the  magnitude  of  is  significant  the  percent  error  is 
relatively  small.  However,  when  the  magnitude  of  Id  approaches  0  the  percent  error  can 
be  quite  large.  The  latter  effect  is  due  to  an  error  in  the  phase  of  the  approximation  with 
respect  to  the  exact  solution.  This  effect  will  occur  when  the  approximation  does  not  have 
the  same  zero  crossings  as  the  exact  integral.  This  applies  to  any  approximation  technique 
and  is  therefore  not  specific  to  the  present  procedure.  Note  that  where  the  approximation 
is  matched  at  a  zero,  i.e.  at  infinite  refractive  index,  the  percent  error  is  small.  Evidently 
matching  one  or  more  zeros  of  the  integral  could  be  a  useful  criterion  in  the  choice  of  an 
adequate  form  of  h.  However,  in  the  present  case  the  quantity  of  physical  interest  is  Qd 
and  thus,  whenever  Id  is  small  the  percent  error  in  Qd  will  be  also  small.  For  example, 
for  the  worst  case  in  Table  I,  m  =  2.0  and  s  =  5,  where  the  error  in  Id  is  103.7%,  the 
corresponding  percent  error  in  Qd  is  only  11.0%. 


The  integral  in  this  example  is  taken  from  (Ref.  6)  and,  as  in  the  previous  example, 
occurs  in  the  scalar  theory  of  electromagnetic  and  acoustic  scattering.  This  integral  accounts 
for  the  effect  of  edge  scattering  from  randomly  oriented  prolate  spheroids  of  arbitrary  aspect 
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TABLE  I 

Error  table  for  Example  1 


m 

■ 

h 

Exact 

Id 

Approximation 

— %  Error — 

1.1 

.001 

1-. 0003129021 

l-.000331086i 

0,5.8 

1.1 

.01 

.999994-. 003130151 

.999993-.00331084i 

.0001,5.8 

1.1 

.1 

.999439-. 03129431 

.999337- .03309771 

.01,5.8 

1.1 

1 

.944813-.30528/?! 

.935227-.3204151 

1. 0,5.0 

1.1 

2 

.790095- .56658ri 

.758557-..581541i 

4.0,2.6 

1.1 

5 

.0551996- .849*2201 

.0305116-.781481i 

44.7,8.0 

1.1 

10 

-.592805- .2373941 

-.540190- .2722221 

8.9,14.7 

1.1 

50 

.0609117-1-.1978221 

.0487599-t-.217617i 

19.9,10.0 

1.1 

100 

-.0911145-.0580755i 

-.0989528- .072956 11 

8.0,25.6 

1.5 

1 

.258516- .9323351 

.226920-.925756i 

12.2,0.7 

1.5 

2 

-.743586- .4655331 

-.719633-.407902i 

3.2,12.4 

1.5 

5 

.490254-1-. 090491 

.416438+. 1620191 

15.1,79.0 

1.5 

10 

.0836259-1- .2696681 

.0418624+. 2760251 

49.9,2.4 

1.5 

50 

.0192536- .05681621 

.0215849-. 06376071 

12.1,12.2 

1.5 

100 

.0160325-.0253829i 

.0183567-. 03078231 

14.5,21.3 

1.5- .Oli 

1 

.25309-.911135i 

.222336-.904551 

12.2,.72 

1.5-.01i 

10 

.065439-f  .2205191 

.0320381+.2262461 

51.0,  2.6 

1.5-.01j 

100 

.00221063-.003401891 

.00259562-. 004020521 

17.4,18.2 

iOl 

1 

.207354-.741218: 

.183119-.734624i 

11.7,.89 

10 

.00677562-1- .035394 11 

.00220001+.03889461 

67.5,4.2 

1.5- li 

1 

.0183699-.1071241 

.0152086- .0145821 

17.2,2.4 

2.0 

1 

-.700659-  6413351 

-.706793- .6037801 

.88,5.9 

.0469974-t-.813015i 

.0734549+. 7433661 

56.3,8.6 

2.0 

.0958382-1- .3456851 

-.00351873+.373174i 

103.7,8.0 

2.0 

10 

-.Ib8222-.10S5861 

-.164505-.162528i 

4.0,48.3 

2.0 

50 

.0216164-.033816ri 

.0250502-. 04264321 

15.9,26.1 

2.0 

100 

.0175966-.00939819i 

.0225043- .009500061 

27.9,1.1 

35-35i 

.001 

.930190-. 06396691 

.930189- .06401051 

.0001,.068 

35-35i 

.01 

.384032- .3147631 

.383882- .3149411 

.04,.057 

35-35i 

.05 

-.0289054-f  .008669231 

-.0288829+. 008736021 

.08, .77 

35-35i 

.09 

.00182362-f  .0001900881 

.00182381+.0001823451 

.01,4.1 
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ratio. 

[19] 

Jo 

where  a,  r  and  co  are  constants  and  P  —  [co8^(ff)  +  r*  8in’(tf)|*/^.  For  this  case  r,  the  aspect 
ratio,  can  vary  from  1  to  oo,  j,  again  the  size  parameter,  can  vary  from  0  to  oo  and  the 
Fock  constant,  cq  »  1. 

Identify  h  with  2Fi[-§,  5;  1;  1  -  5(Ai,A)  with  /with  p8in(d),  a  =  0 

and  b  =  t/2.  For  the  parameters  we  have  Ai  =  cn(r/3)j,  and  771  =  <Ti  =  r.  Thus  we  have 
the  simple  set  p  =  {Ai.r^i}. 

Let 

J(y,)  =  rpsinid)de 
Jo 

=  :[vrr7?t  !ia^]  -  r[co.(,.)v/r:7;;^+ 

and 

=  1  -  l/r^  [20] 

Figure  1  is  a  plot  of  h{9)  and  J{9)  versus  9  for  various  values  of  r.  These  curves 
are  very  similar  and  it  is  evident  that  there  should  be  a  simple  transformation  that  would 
take  one  function  into  the  other 

We  find  d  =  2  a  reasonably  good  approximation  to  the  curvature  of  J  if  we  use  the 


same  functional  form  for  u  as  in  the  previous  example.  As  before,  this  form  gives  the 
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maximum  number  of  three  fitt-'n.  „ 

‘aree  nttmg  parameters,  i.e  n  -  l 

’  ”  “  again  we  find 

J(vi)  Sis  —  ft 

^  +  /=  1,2,3. 

A.  Wor.,  ,h,  ^ 

W  aeep  the  equations  simple  i  e  n 

«dto  =  T/2.  Withthecho«a«.tof«  ^  ‘  -  0,  =  t/4 

VI  and  using  the  fact  that  h(0)  -  t  tj, 

{21]  is  W  -  1,  the  solution  to 


[21] 


ft  = 

ft  = 

ft  = 

7  = 


7ftW!r/2)-q3] 

7ft-A(x/2) 

^A(T/2),hfr/4^ 

7  +  (1  -  «(T/4)J/y(T/47 

A(y/2)  -  I 
•/(»/2)  ' 


(22] 
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TABLE  n 


Error  table  for  Example  2 


■ 

8 

le 

.  Exact 

h 

Approximation 

— %  Error — 

1 

8 

h 

/e 

0 

1.25 

.01 

3.65592x10-® 

3.10578x10-* 

15.05 

1.25 

.1 

.0153844 

.0149671 

2.71 

1.25 

1 

.456008 

.453797 

0.48 

1.25 

2 

.645879 

.643957 

0.30 

1.25 

5 

.847290 

.845950 

0.16 

1.25 

10 

.954779 

.953836 

0.10 

1.25 

100 

1.12007 

1.11984 

0.02 

1.5 

.01 

1.92704x10-® 

1.57466x10-* 

18.29 

1.5 

.1 

.0230733 

.0222162 

3.71 

1.5 

1 

.549222 

.545706 

0.64 

1.5 

.764053 

.761085 

0.39 

1.5 

.989166 

.987146 

0.20 

1.5 

1.10846 

1.10705 

0.13 

1.5 

1.29101 

1.29067 

0.03 

2.0 

1.45147x10"^ 

1.20570x10"^ 

16.9 

2.0 

n 

.0419768 

.0404307 

3.68 

2.0 

1 

.745312 

.740096 

0.70 

2.0 

1.01026 

1.00589 

0.43 

2.0 

1.28302 

1.28007 

0.23 

2.0 

10 

1.42612 

1.42407 

0.14 

2.0 

100 

1.64362 

1.64311 

0.03 

5.0 

.01 

3.89743x10-* 

3.61508x10"* 

7.24 

5.0 

.1 

.385611 

.380045 

1.44 

5.0 

1 

1.97382 

1.95165 

1.12 

5.0 

2.55329 

2.53113 

0.87 

5.0 

3.12749 

3.10992 

0.56 

5.0 

10 

3.42199 

3.40885 

0.33 

5.0 

100 

3.8627 

3.85911 

0.09 
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TABLE  II  (continued) 


r 

8 

/. 

Exact 

/e 

Approximation 

— %  Error — 

10 

.01 

1.22907x10-* 

1.17944x10-® 

4.04 

10 

.1 

.385611 

J80045 

1.44 

10 

1 

4.02320 

3.96007 

1.57 

10 

2 

5.14791 

5.08044 

1.31 

10 

5 

6.24900 

6.19179 

0.92 

10 

10 

6.80926 

6.76488 

0.65 

10 

100 

7.64262 

7.62982 

0.17 

100 

.01 

.00014741 

.000144273 

2.12 

100 

.1 

4.03452 

3.97989 

1.35 

100 

1 

40.6192 

39.7725 

2.08 

100 

2 

51.7342 

50.7669 

1.87 

100 

5 

62.5487 

61.6353 

1.46 

100 

10 

68.0260 

67.2530 

1.14 

100 

100 

76.1406 

75.8643 

0.36 

lO-* 

.01 

.0147756 

.014492 

1.92 

lO** 

1 

4062.54 

3974.18 

2.17 

10* 

100 

7613.79 

7576.89 

0.48 

10* 

.01 

.147756 

.144931 

1.91 

10* 

100 

76137.9 

75764.8 

0.49 

10* 

.01 

1.47756 

1.44923 

1.92 

10* 

100 

761379. 

757639. 

0.49 

With  the  above  equations  we  can  substitute  into  [4]  to  get: 

fff  g-Xik 

=  .i.k 

a  s  h(0)  =  1 


[23] 


0  s  /»(x/2). 
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Tabh  n  presents  the  error  between  the  exact  numeric;Jly  calcnlated  value  of  and  [24]  for 
a  wide  range  of  r  and  s  values. 

For  an  aspect  ratio  of  1,  the  approximation  is  equal  to  It  exactly  for  any  size 
parameter  s.  The  error  is  greatest  for  a  given  r  when  s  is  small.  This  is  once  again  due  to 
the  presence  of  a  zero  in  the  integral  for  s  =  0.  The  approximation  is  exact  (i.e.,  equal  to  0) 
for  s  =  0,  however  the  asymptotic  behavior  as  s  approaches  0  is  not  quite  correct.  This  is 
not  important  for  the  physical  problem  of  edge  scattering  since  the  real  quantity  of  interest 
is  related  to  2  -  /e(s)//t(oo). 

3.3  Example  3 

This  example  occurred  in  the  study  of  scattering  from  randomly  oriented  infinite 
cylinders  (Ref.  7). 

r^^Uiiw)sia\0)d9  [24] 

Jo 

where  Hi  is  the  first  order  Struve  function  and  to  =  2j[(nj*  -  cos^(tf))*^*  -  8in(0)].  The 
size  parameter  a  can  vary  from  0  to  oo,  Re(m)  can  vary  from  1  to  oo,  and  Im(m)  from  0  to 
'•00.  Again  m  is  the  index  of  refraction. 

In  this  integral  h  =  w/2a,  g(Xi,h)  =  Hi(A|/i),  f{9)  =  sin’(tf),  o  =  0,  and  b  =  t/2. 
The  parameters  are  Ai  =  2j  and  r/i  =  m.  Thus  p  =  {A],f/i}.  Again  the  integrabiiity  con¬ 
dition  [4]  must  be  satisfied.  Since  s(Ai/ij  =  Hi(AiA)  =  iFj[l;3/2,5/2; -Aj/»^/4],  we  can 
transform  [24]  into  an  integral  of  a  form  similar  to  [9].  If  we  consider  [10],  the  hypergeomet- 
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ric  representation  of  Hi  as  gven  above  and  in  [9],  then  the  choices  for  u  will  be  “restricted” 
to  a  general  hypergeometric  pF^.  However,  since  attention  must  be  given  to  the  compu¬ 
tational  feasibility  of  the  solution,  only  the  lowest  order  hypergeometric  functions  should 
be  considered.  Note  also  that  the  exponent  of  the  variable  h  in  the  chosen  hypergeometric 
function  must  be  2  in  accordance  with  [9]. 

Let 

J(yi)  = 


/■ 


sin*(d)  d9 


pj  _  cos(yi)5in(yi) 
2  2 


[25] 


Figure  2  is  a  plot  of  h(9)  and  J(9)  versus  9.  Again  as  in  the  previous  example, 
these  curves  are  similar  in  that  they  are  both  monotonic  with  gentle  curvature.  Given  the 

t 

restrictions  placed  on  u  as  detailed  above,  the  simplest  hypergeometric  function  that  can 
be  used  is  o.fb(;cA^]  =  In  this  case  the  approximation  to  J  would  be 


J'(yj)  »  ?iErf(93A(y,)]  +  ft,  /  =  1,2,3. 


[26] 


Choosing  yi  =  0,  yj  =  1  and  ys  =  t/2  and  solving  for  the  fitting  parameters  we 

get: 


_ _ 

Erf{9jA(0)]  -  Erf[ftA(T/2)] 


ft  =  “9i  Erf(ftA(0)] 


[27] 
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FIGURE  2  -  A  plot  of  h($)  ajid  J{6)  versus  0  (Example  3) 
aod  93  is  the  root  of 


Erfl^aMO)]-  Erft93Ml))  _  1 ,, ,  x 
Erfl92/»(0))  -  Erf{g2M»/2)]  »  ^  ^  ' 

Here  we  have  the  situation  where  one  of  the  fitting  parameters,  in  this  case  93,  cannot  be 

solved  for  explicitly.  We  therefore  must  approximate  the  solution  of  [281  with  some  simple 

analytic  function.  This  can  be  done  by  replacing  the  Erf  functions  in  [28]  v/ith  their  first 

order  asymptotic  expansion.  This  is  justified  since  it  can  be  shown  that  qih{yi)  is  always 

greater  than  1.49.  Doing  this,  we  find  that  the  limiting  function,  both  for  large  and  small 

Abs(m  -  1)  is: 


1  r4j(i)h(i)' 


[29] 


Since  the  asymptotic  expansion  has  the  worst  error  for  small  Abs(m  -  1),  all  other  values 
of  m  will  produce  a  better  approximation.  When  the  above  approximation  is  used  an  error 
of  1%  or  less  is  introduced  In  the  calculation. 


i 

t 
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Table  HI  is  a  list  of  the  fitting  parameters  q^  and  q2  with  h{Q)  and  h(5r/2)  for 


vanous  m. 


Table  of  fitting  parameters 


h(T/2) 


1.3-i 

1.5 

1.5- .01i 

1.5- .li 
1,5-i 


10.0 

100.0 

1000.0 


•27.9301 

-34.0891 

-74.1495 

-74.0973-1-2.833371 
74.1207-46.9499  i 
•  157.782 

-157.676+5.89236i 

-147.167-t-57.6607i 

155.283-89.37981 

-987.166 

-34496.5 

-1.12303x10'* 

-3.52471x10^ 

-1.08171x10* 

-3.26993x10'° 

-9.80490x10" 

-3.10274x10'° 

-8.53632x10'^ 

-3.5921x10'" 

•  1.80814x10'"^ 


31.0506 

16.07096 

6.02426 

6.0187-i-.169276i 
1.32201-fl.60465  i 
3.97178 

3.97057-f  .06240251 
3.85538+. 6011661 
1.48387+1.398531 
2.37186 
1.4976 
1.16944 
.9887.50 
.871188 
.787146 
.723318 
.672757 
.631446 
.184026 
.05n47 


.320156 

.458257 

.83066 

.830750- .01564851 
1.07434-1.21004  i 
1.118035 
1.11807..013416i 
1.12156-.133743i 
1.27679-1.174821 
1.73205 
2.82843 
3.87298 
4.89898 
5.9161 
6.92820 
7.93725 
8.94427 
9.94987 
99.995 
999.999 
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Applying  the  method  to  [24]  with  the  above  information  we  have: 


/»  ss  - . ^  f 

V*’  A(0) 

~  3x3/2  A(0)  ^  ‘^V/2,  5/2  4  j 


491^2  A? 


3/2 :  1; 

5/2:  3/2,  5/2; 


2  r 

«  l»/<Jlfi:2;o|^5/2:  3/2,  5/2; 


,,, 


,-qlh^{rl2) 


]} 


[30] 


where  in  the  last  equation  the  Kampe  de  Feriet  functions  have  been  reduced  from  Fjijio  to 
^!:2;0  since  a  set  of  coefficients  coincided.  Simple  computer  routines  for  evaluating  Kampe 
de  Feriet  functions  are  given  in  (Ref.  4).  These  were  used  in  slightly  modified  form  to 
compute  Table  IV  for  small  values  of  s.  Other  more  efficient  algorithms  for  computing  [30] 
are  given  in  Appendix  A.  These  algorithms  are  derived  from  the  concepts  found  in  Refs.  8 
and  9. 


Table  FV  gives  errors  that  are  generally  less  than  10%  for  real  m  since  the  only  zero 
of  the  integral  occurs  at  s  =  0  as  in  Example  1.  For  complex  values  of  m  and  sufficiently  large 
values  of  s,  I]  becomes  an  exponentially  increasing  function  of  s.  So  does  the  approximation 
[30].  For  any  two  rapir'ly  increasing  functions,  a  slight  displace  -ent  in  the  independent 
variable  of  one  of  the  functions  will  produce  a  substantial  percentage  difference  between  the 
values  of  the  two  functions.  Since  [30]  has  a  slight  difference  in  phase  (i.e.,  displacement  in 
j)  in  comparison  with  /‘,  large  percent  errors  occur.  This  can  be  seen  in  the  error  table. 


UNCLASSIFIED 

22 


TABLE  IV 

Error  table  for  Example  3a 


m 

s 

n 

Exact 

n 

Approximation 

— %  Error — 

1.1 

.0025896 

.00223957 

13.5 

1.1 

.0103105 

.00893241 

13.4 

1.1 

.223644 

.203632 

8.9 

1.1 

10 

.608131 

.608241 

0.02 

1.1 

20 

.708387 

.722618 

2.0 

1.1 

50 

.595284 

.573136 

3.7 

l.l-.li 

.5 

.00044929- .005016631 

.000147684-. 004464611 

67.1,11.0 

l.l-.ii 

1 

.00194818-.0200093i 

.000692555- .01785041 

255.5,10.8 

l.l-.li 

.151968-.444002i 

.0974261-.432831i 

35.8,2.5 

1.3 

.0204227 

.0190081 

6.9 

1.3 

1 

.0793354 

.0742178 

6.5 

1.3 

5 

.808644 

.832955 

3.0 

1.3 

.359062 

.323648 

9.9 

1.3 

20 

,431977 

.410225 

5.0 

1.3 

.496680 

.482361 

2.9 

1.3- .Oli 

.5 

.0204049-. 00126728i 

.0189901-.00121746i 

6.9,3.9 

1.3-.01i 

1 

.0792763-. 00478980i 

.0741567-.00463915i 

6.5,3. 1 

1.3-i 

.5 

-.156730-.157030i 

-.162314-.153220i 

3.5, 2.4 

1.3-i 

1 

-.610466- .8612051 

-.652805- .8554131 

6.9,0.7 

1.3-i 

1244.224976-I-951 .90549 17i 

1588.029043-1-1644.5693081 

27.6,72.8 

1.5 

.0525496 

.0502753 

4.4 

1,5 

1 

.195683 

.188589 

3.6 

1.5 

5 

.512050 

.503202 

1.7 

1.5 

.611504 

.626669 

2.5 

1.5 

il 

.471594 

.499574 

5.9 

1.5 

il 

.595284 

.573136 

3.7 

1.5-.01  i 

.5 

.0525342-.00192994  i 

.0502600-. 00189406  1 

4.3,1.9 

1.5-.01  i 

1 

.195648-. 00668855  i 

.188554- .00664175  1 

3.6, 0.7 

1.5-.01  i 

.613522-I-.011374  i 

.611768-f  .0187999  1 

0.3,65.2 

1.5-.1  i 

,5 

.0510075-. 01933631 

.0486982-. 01897801 

4.5,  1.8 

1.5-.1  i 

1 

.192127-. 0672897  i 

.184783-. 0668288  1 

3.8,0.7 

1.5-.1  i 

5 

.4737- .2601  i 

.4571-. 2974  1 

3.5,  14 

1.5-i 

.5 

-.106388-. 227092  i 

-.122096- .224601  i 

14.8,1.1 

1.5-i 

1 

-.218598-1.12719  i 

-.254697-1.13575  1 

16.5,0.8 
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I 


TABLE  rV  (continued) 


m 

s 

n 

Exact 

n 

Approximation 

— %  Error — 

.182774 

.179493 

1.8 

Bl 

1 

.567054 

.562553 

2.5 

2.0 

5 

.629007 

.640619 

1.8 

Bl 

.466074 

.488837 

4.9 

m 

.522300 

.527632 

1.0 

Bl 

.488015 

.492513 

0.9 

.5 

.181839- .03182501 

.178535-.0317615i 

1.8, 0.2 

1 

.569583-. 07178921 

.565011-.0730041i 

0.8, 1.7 

.5 

.822542 

.824195 

0.2 

5.0 

1 

.434339 

.470611 

8.4 

5 

.53203 

.547058 

2.8 

5.0 

10 

.50478 

.507859 

0.6 

20 

.513551 

.490891 

4.4 

50 

.502995 

.500390 

0.5 

An  additional  advantage  of  this  technique  is  that,  once  the  work  has  been  done  to 


find  a  satisfactory  approximation  such  as  [2]  and  [3],  other  integrals  are  readily  solved  with 
different  g{u)  to  a  similar  level  of  accuracy.  For  example,  if  the  Hj  function  of  the  above 


example  is  changed  to  a  first  order  Bessel  function  Ji,  the  integral  [24]  becomes: 


/’=  r^^7i(AiA)  sin*(tf)d<?. 

Jq 

With  the  same  substitutions  as  used  to  obtain  the  solution  to  [24]  we  get: 

A(o) 


[31] 
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=  -  W.(^,A.MO))] 

))  -  (32] 

where  U\  and  U3  sire  Lommel  functions  of  two  variables  (Refs.  8  and  10). 

Table  V  is  the  error  table  for  this  approximation.  Unlike  with  Hi,  this  integral 
has  zeros  for  real  values  of  m  and  s  ^  0.  Thus,  the  phase  problem  as  previously  mentioned 
in  other  examples  occurs.  Again  numerical  schemes  for  evaluation  of  [32]  are  detailed  in 
Appendix  A.  For  complex  m  and  large  s,  the  large  percent  errors  occur  for  the  same  reason 
as  already  explained  for  Table  IV. 

The  quantity  of  physical  interest  in  cylindrical  scattering  theory  is  Qci  the  extinction 
efficiency.  In  terms  of  the  above  symbols,  we  have 

Qc  =  4Re  (/"  +  ill)  ■  [33] 

This  function  has  no  zeros  except  for  3  =  0.  Table  VI  is  an  error  table  of  Qc  for 
complex  m.  For  real  m,  the  percent  error  is  identical  to  that  of  Table  IV.  It  is  now  clearly 
seen  that  the  errors  for  complex  m  are  of  the  same  order  as  for  real  m  even  for  large  3. 
Furthermore,  Qe  approaches  2  for  large  3  and  any  m,  as  predicted  by  Babinet’s  principle  of 
diffraction  from  large  objects. 


I 
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TABLfc  V 

Error  table  for  Example  3b 


m 

— 

s 

Exact 

n 

Approximation 

— %  Error — 

1.1 

.5 

.0472756 

.0451376 

4.5 

1.1 

1 

.0939029 

.0898062 

4.4 

1.1 

5 

.377706 

.378477 

0.2 

1.1 

.381027 

.413833 

8.6 

1.1 

20 

-.126605 

-.166460 

31.5 

1.1 

-.0425576 

-4)78062 

83.4 

l.l-.li 

.5 

.04858 12- .04534.'i5i 

.0459329-.0443288i 

5.5, 2.2 

l.l-.li 

1 

.0981060-.0S92763i 

.0927459- .08767081 

5.5, 1.8 

l.l-.li 

.591899-. 217036i 

.583436-.265611i 

1.4,22.4 

1.3 

.133658 

.130523 

2.3 

1.3 

1 

.254182 

.249666 

0.5 

1.3 

5 

.136981 

.150056 

9.5 

1.3 

10 

-.6840142 

-.0928197 

10.5 

1.3 

20 

-.0591054 

-.0296963 

50.2 

1.3 

50 

-.0376139 

-.0206433 

54.9 

1.3- .Oli 

.5 

.133668-.00410752i 

.130533-.00409969i 

2.3,0.2 

1.3-.01i 

1 

.254235-. 00702857i 

.249718-.007119171 

1.8,1.3 

1.3-i 

.5 

,219293-.441338  i 

.213943-.446550i 

2.4,1.1 

1.3-i 

1 

.912076-1.0172801 

.905500-1.058881 

0.7,4.1 

1.3-i 

5 

-951.902855-f- 1243.7280891 

-1644.566682-f-1587.532198i 

72.8,27.6 

1.5 

.5 

.211668 

4»8736 

0.3 

1.5 

1 

.371618 

.369804 

0.5 

1.5 

5 

-.215697 

-.241939 

12.2 

1.5 

20 

.0528324 

4)570300 

7.9 

1.5 

50 

.00725979 

.00539024 

74.2  1 
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TABLE  V  (continued) 


m 

B 

Exact 

n 

Approximation 

— %  Error — 

1.5- .Oli 

.5 

.211679-.00368736i 

.208748-. 003710241 

1.4,0.62 

l.S-.Oli 

1 

.371684- .00463889i 

.369871-.00480011i 

0.48,3.5 

1.3- .Oli 

10 

-.04818764-.0211236i 

-.0831952+ .02055741 

72.6,2.8 

1.5-.li 

.212773-. 03689711 

.209800- .03713271 

1. 5,0.6 

1.5-.li 

.378135- .04653031 

.37627-.0481714i 

G.5,3.5 

1.5-.li 

-.330915- .04913881 

-.376455- .06287561 

13.8,27.9 

1.5-i 

.320290- .402570  1 

.316487-. 407628  1 

1.2,1.3 

1.5-i 

1 

1.19943-.643665  1 

1.20758- .678831  1 

0.7,5.5 

mm 

.5 

.365144 

.363883 

0.3 

1 

.433180 

.438278 

1.2 

5 

-.0496403 

-.0725961 

46.2 

2.0 

El 

.0602173 

.0729676 

21.2 

ml 

.0274814 

.00953082 

65.3 

2.0 

-.00742132 

.00106124 

114.3 

.5 

.366629-. 023893ri 

.365371-.0242659i 

0.3, 1.6 

1 

.439381-|-.0214225i 

.444627+. 0205887 

1.2,3.9 

-.0866975 

-.114949 

32.6 

.197090 

.206872 

5.0 

.0452478 

.0417831 

7.7 

10 

-.0275163 

-.0251509 

8.6 

5.0 

20 

.00346559 

-.00354774 

202.4 

in 

50 

-.00472154 

-.00248860 

52.7 
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TABLE  VI 


Error  table  of  Qc 


m 

■ 

Exact 

Approximation 

— %  Error — 

l.l-.li 

.5 

.1831792 

.1779059 

2.9 

l.l-.li 

1 

.3648979 

.3534534 

3.1 

l.l-.li 

1.476016 

1.452150 

1.6 

1.3- .Oli 

.0980497 

.0923592 

5.8 

1.3-.01i 

1 

.345219 

.325103 

5.8 

1.3-i 

.5 

1.138432 

1.136944 

0.1 

1.3-1 

1 

1.627256 

1.624300 

0.2 

1.3-i 

1.987548 

1.987380 

.008 

1.5- .Oli 

.224886 

.215881 

4.0 

1.5-.01i 

1 

.801148 

.773416 

3.5 

1.5-.01i 

10 

2.369586 

2.364842 

0.2 

1.5-.li 

.5 

.351618 

.343325 

2.4 

1.5-.li 

1 

.954629 

.931818 

2.4 

1.5-.li 

2.09136 

2.07990 

0.5 

1.5-i 

1.184728 

1.142128 

3.6 

1.5-i 

1 

1.700268 

1.696536 

0.2 

»«2i| 

.5 

.8229108 

.8112036 

1.4 

1 

2.192642 

2.177689 

0.7 

4.0  INTEGRALS  INVOLVING  THE  EXPONENTUL  FUNCTION 


In  this  chapter,  we  will  indicate  how  this  integration  technique  allows  for  a  trade-off 
between  the  complexity  of  the  special  function  used  and  the  accuracy  of  the  solution.  Since 
several  of  the  examples  shown  in  the  last  chapter  involve  integrals  owr  the  exponential 
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function  of  a  complicated  argument,  we  will  use  this  form  of  integral  to  demonstrate  this 
point.  (It  is  not  coincidental  that  integrals  over  the  exponential  function  are  common  since 
they  are  specific  caa<»«  cf  general  diffraction-type  integrals.) 


Consider  the  integral 


The  technique  then  requires 


[35] 


If  one  chooses  the  form  of  u{h)  to  be  a  polynomial  of  h  with  integer  powers,  the  resulting 
approximation  to  I  will  be  in  terms  of  a  sum  of  exponential  integrals.  Thus 


where  m  and  n  are  any  positive  integers  or  zero.  If  we  chose  u(h)  to  be  a  polynomial 
of  h  with  rational,  real  or  complex  powers,  then  the  approximation  becomes  a  series  of 


incomplete  gamma  functions.  Hence, 


=  f;  Ok  (r[/i(ib)  -b  l,/»(o)|  -  rlM(k)  +  1,  A(6)])  h{a)  and  h{b)  ^  0 

n 

=  H  (tW*)  +  l.M*)]  -  7[/‘W  +  1.  A(®)])  R«[m(^)]  +  1  >  0.  [37] 

iatO 


Obviously,  as  one  increases  the  number  of  terms  in  the  polynomial  representation  of 
u(h),  the  accuracy  of  the  solution  will  be  improved  at  the  cost  of  adding  additional  terms. 
This  allows  for  balancing  solution  complexity  against  accuracy. 
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5.0  CONCLUSIONS  AND  DISCUSSION 

We  have  developed  a  method  of  obtaining  an  analytic  appioximation  to  integrals 
that  may  have  complicated  functions  in  their  kernels.  The  demonstrated  advantages  of  the 
technique  are: 

•  sometimes  producing  a  compact,  easily  calculable  form; 

•  approximations  to  integrals  where  standard  approaches  fail; 

•  possible  trade-offs  between  functional  complexity  of  the  solution  against  accuracy; 
and 

•  solutions  that  are  often  globally  valid  over  the  full  range  of  each  parameter. 

Not  all  of  the  above  advantages  may  occur  when  using  the  technique.  When  compact  forms 
can  be  achieved  not  only  can  the  integral  then  be  rapidly  computed  but  considerable  gmns 
in  physical  insight  can  follow.  Having  the  ability  to  compute  rapidly,  even  approximately, 
difficult  integrals  is  a  fundamental  requirement  in  design  and  inverse  problems  where  pa¬ 
rameters  have  to  be  varied  over  large  ranges  and  statistics  established  for  the  sensitivity  of 
the  results. 

We  are  using  this  integration  technique  to  solve  many  integrals  involved  in  the 
evaluation  of  the  extinction  efficiency  of  randomly  oriented  discs,  infinite  cylinders  and 
prolate  and  oblate  spheroids  for  arbitrary  materials  and  sizes. 
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Since  the  form  of  these  integrals  is  the  same  as  those  of  general  diffraction  theory,  we 
expect  the  technique  to  be  applicable  to  many  problems  in  optics,  as  for  example  the  design 
of  Raman  cells  using  exotic  laser  pumping  geometries  such  as  axicons.  We  also  expect  the 
technique  to  be  useful  in  other  fields.  It  is  at  this  stage  a  possible  solution  looking  for  other 
promising  problems. 
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APPENPK  A 

Namerical  Schemes  for  Evalnatiag  Examples  3a  and  3b 


Evaluating  the  Feriet  functions  in  [30]  and  [32]  directly  as  double  sums  works  only  if  both 
the  arguments  are  small.  K  one  or  both  are  large,  the  series  will  still  converge  but  the 
intermediate  sums  will  be  extremely  large.  When  this  occurs,  multiprecision  calculations 
would  be  required  for  the  converged  answer  to  retain  any  significance.  At  this  point  the 
calculation  becomes  impractical.  Hence  other  approaches  are  necessary. 


An  algorithm  that  works  for  small  values  of  Xih  can  be  derived  from  the  series 
expansion  for  and  Jk(z)  (Ref.  8): 


H„(r)  + 


OP 


E 


r(m/2+l)r(rj/2+H.i^)- 


[A  1] 


Performing  the  integration  as  indicated  by  [30]  and  [32]  on  [A  1]  term  by  term  with  i/  =  1 
gives: 


ri,-r3  _  f  h(T/2)  ^  e»0-"»W»[A,A(x/2)/2r-H7-(m/2  +  1,  g|/i»(T/2)) 

+  ~  \  2  r(m/2  +  2) 

2  ^  r(m/2  +  2)  / 

Evidently  m  odd,  [A  2]  ^ves  the  expansion  for  /'  and  for  m  even,  the  series  for  It 
is  clear  that  large  Xih  will  create  the  same  summation  difficulties  as  it  does  for  the  direct 
calculation  of  the  Feriet  functions.  For  the  case  of  large  Ajh,  we  can  derive  an  asymptotic 
series  for  by  a  similar  term  by  term  integration  as  was  done  for  [A  2]  but  over  the 
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asymptotic  series  for  (Ref.  8).  Thus  taking  the  asymptotic  series 

H  I  ^  V  r(k  +  i/2) 

H^(^)  -  Yu{2)  +  ,  ^  r(j/  +  1/2  -  Ifc)(f 

we  get,  again  setting  v  =  1, 


[A  3] 


rl  1  ,  2^193 


h(0)  ^  r(Jb  +  i/2)Ek+i/2[gMo)] 
T  ^  r(3/2  -  k)[Ajh2(0)]* 


k=l 


hiir/2)  A  Tik  +  l/2)E,^^,,[qlhHrf2)]  ,  291ft  /*<*/2) 


*  ■A  A(oi 


[A  4] 


T  T(Z/2-k)[Xih\ir/2)]^  '  Jk(o) 

The  integral  in  [A  4]  is  the  same  as  [32]  with  Vi  replaced  by  J\.  Since  the  properties  of 
Vi  and  Ji  are  very  similar,  both  integrals  can  be  evaluated  with  a  procedure  by  Agrest 
and  Maksimov  (Ref.  9,  p.l23ff).  They  develop  two  series  for  either  or  FJ,  according  to 
whether  2q\hlX\  <  1  or  >  1.  In  both  cases,  multiple  integration  by  parts  is  performed 
using  basic  properties  of  these  functions.  The  details  can  be  found  in  the  cited  reference. 
The  results  are  for  J\: 

J2  «  2ftft  f 


29^(0) 


A[AiA(0)] 


for  293 A/ Ai  <  1  and 


ja  „  29193  I  . 


k>0 


l*+l 


29|/>(t/2) 

Ai 


A[A,A(t/2)]  [A  5] 


for  2q\h/X\  >  1, 
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and  letting  stand  for  the  analogous  case  for  the  integral  over  Y\,  we  get: 


for  2q\hlX\  <  1  and 

oo  r  \  1*+1  I 

for  2qlhjXi  >  1. 


An  additional  numerical  scheme  for  computing  [32]  by  using  the  double  variable 


Lommel  function  representation  can  be  derived  from  the  appendix  of  Ref.  10. 
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